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Abstract 

A Monte Carlo simulation is carried out to investigate the magnetic ordering in magnets with 
random anisotropy (RA). Our results show peculiar similarities to recent experiments that the real 
part of ac susceptibility presents two peaks for weak RA and only one for strong RA regardless 
of glassy critical dynamics manifested for them. We demonstrate that the thermodynamic nature 
of the low-temperature peak is a ferromagneticdike dynamic phase transition to quasidong range 
order (QLRO) for the former. Our simulation, therefore, is able to be incorporated with the 
experiments to help clarify the existence of the QLRO theoretically predicted so far. 

Keywords: Dynamic transition, random magnetic anisotropy, amorphous magnets, Monte Carlo simulation, 
critical slowing down, complex susceptibility 



*Author to whom correspondence should be addressed. Electronic mail: nmha@ess.nthu.edu.tw 
^Electronic mail: pyhsiao@ess.nthu.edu.tw 



1 



Over decades, great interest has been intensively addressed to rare-earth magnetic glasses 
of random magnetic anisotropy (RMA). However, the nature of magnetic phase transition 
(MPT) and magnetic ordering in such random magnets has still been far from being com- 
pletely understood [1]. First, Monte Carlo (MC) simulations were opposite to the prediction 
of renormalization group theories to reveal that the second-order MPT exists in three- 
dimensional weak RMA systems of XY [4] and Heisenberg spins {5], (|. Second, Itakura and 
Arakawa 9] have demonstrated that a crucial additional vortex energy should be included 
in the Imry-Ma type arguments, which have predicted the absence of ferromagnetic long 
range order (LRO) in magnets of random field (RF) [7] and RMA 8] for space dimensions 
d < 4, to explain the power-law correlation of quasi-long range order (QLRO) in the Bragg 
glass state of impure superconductors [10], and showed MC results of the power-law scenario 



for the weak RF model of XY spins. Feldman [TJJ has theoretically shown that QLRO can 
emerge instead of LRO in d = 4 — e dimensions and is common in such impure systems of 
continuous non-Abelian symmetry as magnets of weak RF and RMA. In addition, QLRO 
has also been clearly evidenced in a number of MC simulations [a, |9|, 1121], of which a power- 
law spin correlation function has been found to indicate a ground state of QLRO in weak 
RMA systems of Heisenberg spins [5]. In spite of these theoretical conjectures, the lack of 
direct experimental and theoretical agreements in the literature leads to the questions of (i) 
whether the magnetic transition and low-temperature magnetic order in weak RMA mag- 
nets are ferromagnetic-like, and (ii) whether the so-called RMA model [l3] can be applied 
to understand such phenomena in real materials. 

In this Letter, wc address these questions by conducting a MC simulation upon the RMA 
model l^l . Our simulation aims to clarify an experimental possibility that the singularity on 
the temperature-dependent curves of the real part of the ac susceptib ility , x'(T, u), for a weak 
RMA glass of a-Ho2sFe72 amorphous film reported by Saito et al. 14j manifests a second- 
order MPT, which is discriminated in nature from the magnetic glassy phase transition 
(GPT) in strong RMA glasses, for instance, Dy4oAl 2 4Co2oYiiZr 5 bulk glass [l^]. Differing 
from the MPT, the GPT is indicated by a glassy critical slowing down law without any 
singularity shown on x'(T,u) at the transition temperature, T g [15(. Notice that a similar 
scenario, but for a case of RF systems, has been existed in the literature when Schremmer 
and Kleemann [16] demonstrated for an orientational glass system of K4_2.Li3.TaO3 with 
x = 0.063 (a doping well above the glassy and ferroelectric boundary x c ~ 0.022) that the 



2 



singularity on the real part of the ac dielectric permittivity, e'(T, to), is of a transition to the 
long-range ferroelectric phase. However, this is a first-order transition. 

In the light of these experiments, we show in the present work that the nature of the 
aforementioned singularity for the weak RMA magnet of a-Ho2sFe72 amorphous film [14] can 
be understood dynamically with the concept of dynamic transition within the framework of 
the RMA model of three-dimensional Heisenberg spins [3], of which the Hamiltonian can 
be written as 

where the first term is due to the exchange coupling with strength J > between nearest- 
neighbor spins, the second is for on-site RMA with strength D > 0, and the last is the 
Zeeman term with the presence of an external field of strength H along z axis. Si and hi 
are unit vectors representing the spin (an annealed variable) and the random easy axis (a 
randomly- quenched variable) at site i, respectively. In Eq. (pQ), the anisotropy to exchange 
ratio, D/J, plays the role of the degree of RMA. One prominent effect of the degree of RMA 
is clearly observed in Fig. [TJ Her e, V iT, u) is simulated using the same MC technique as that 



described in our previous papers [17|,ll8| for simple-cubic-lattice systems of LxLxL (L = 10) 



Heisenberg spins as an external ac field, H = H sm(ut), is applied, where H /J = 0.05, 
time t is in MC step (MCS), and frequency u is in MCS -1 . Each data point is averaged over 
50 realizations of {hi}f 3 . As shown in Fig. [TJ curves of x'(T,uj) with u = 3 x 10~ 3 exhibit 
two peaks for small values of D/J, i.e., weak RMA {D/J < 5). The high-temperature 
peak is responsible for an Arrhenius-type relaxation which is in common with that for those 
curves of large values of D / J of strong RMA, whereas the low-temperature peak peculiarly 
characterizes another magnetic nature of magnets of weak RMA, the position of which is 
almost insensitive to the change of anisotropy strength. Notice that in our simulation, we 
mimic the measurement protocol that the system is cooled in the ac field to the lowest 
temperature then carrying out the calculation of ac susceptibility and other quantities when 
heating the system up. The reason for this choice is because we have seen in our simulation 
that, unlike the RF system of ¥^i_ x Li x Ta.O-i with x = 0.063 [it] ], cooling the systems of weak 
RMA in a nonzero dc field even as small as Ho shall unexpectedly result in the suppression 
of the low-temperature peak of x'(T,u) curve and the curve looks like that of strong RMA, 
i.e., an one-peak curve. Interestingly, these distinct characteristics of x'(T,ui) for weak and 
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strong RMA systems are consistent with results reported by Itakura that the function 
G(r) oc r^" -1 exp(— r/£) can be used to describe spin correlation of the ground state for the 
RMA model in Eq. (pQ). The correlation length £ is finite for large values of D/J while it 
is infinite for weak RMA of D / J < 5 so that the spin correlation reduces to a frozen power 
law of QLRO ground state, G(r) oc r^ 11 ^ 1 . We remark that we shall only focus on a weak 
RMA glass of D/J = 3.5 and a strong RMA glass of D/J = 10 which are typical of the 



RMA mode 
a-Ho2sFe 7 2 



of Heisenberg spins in Eq. ([I]) to understand magnetic behaviors of weak RMA 



14 1 and strong RMA Dy 40 Al24Co2oY 11 Zr 5 glasses [15|, respectively. 



Figure [2] presents the temperature dependence of x'(T,ui) at different frequencies and 
H /J = 0.05 for D/J = 3.5 and D/J= 10 cases. For D/J = 3.5, curves of x'{T,uj) exhibit 
two peaks. The position of the low-temperature one, T p , is insensitive to frequency and is 
at about T p /J « 1.15. The position of the high-temperature one, T^(uu), shifts toward low 
temperature in an Arrhenius way with decreasing frequency in addition to increasing the 
heights of the two peaks. At sufficient low frequencies, in this case u < 5 x 10" 4 , the two 
peaks merge together so that x'(T,u) rockets up then drops abruptly at about T c / J = 1.0. 
This fashion is what has been experimentally shown for the a-Ho2sFe72 glass and the "dip" 
in x'(T,u), which occurs at the same position of the single peak in x"(T, a;) (not shown), 
indicates a signature of the MPT singularity [14]. Besides, the temperature dependence of 
the ac susceptibility obtained in our MC simulations for the D/J = 3.5 case also shows 
another feature resembling the experiment of a-Ho2sFe72 glass. In contrast to spin glasses 

SGs), x"( u ) > x'iyi) m the vicinity of the singularity, which, according to Saito et al. 

14j , "implies that the center r c of distribution of relaxation time g(lnr) is much longer 
than the measuring time constant t = 1/u." Focussing on the dynamic behavior at the 
transition region, the authors applied a phenomenological Cole-Cole model of polydispersive 
relaxation which yields the ac susceptibility as x( u ) — Xa + (xo — Xa)/{^1 + {i UJT cY} and 
g(lnr) = sin( / 97r)/27r{cosh[/?ln(r/T c )] + cos(/37r)}, where xo an d Xa are static and high 
frequency limit susceptibilities, and < (3 < 1. They found that g(lnr) is almost Gaussian 
in lnr and symmetric about lnr c , (3 reduces from 1 to 0.4 and r c becomes longer and longer 
with decreasing temperature toward T c in company with broadening of g(lnr). All of these 
features are similar to SGs, however, r c for the a-Ho2sFe7 2 glass is several orders of magnitude 
longer than those of SGs. On the other hand, the low-temperature peak is suppressed for 
all frequencies in the D / J = 10 case like that of strong RMA of Dy4oAl 2 4Co2oYnZr5 glasses 
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Q, Q. Nonetheless, we did find that there is a well-determined transition temperature, T g , 
of the GPT for both D / J = 3.5 and D/J = 10 cases by means of the scaling law of critical 
slowing-down dynamics, r c = T*[T^,(uj))/T g — 1]~ ZL/ , shown in the insets of Fig. [21 In terms 
of this scaling law, the magnetic glassy behaviors for systems of weak and strong RMA are 
expected to be the same and like those of SGs 111]. For instance, if the critical exponent 

n n 

v of the correlation length roughly takes values in the range of 0.7 ~ 0.8 |4J, 16[ then the 
dynamical exponent z may be 1.85 ~ 2.35, i.e., consistent with the magnitude of those for 
SGs [19]. Another example is that Billoni et al. [^(J have reported aging phenomena for 
the D/J = 3.5 case similar to those of Heisenberg SGs at low temperatures. To this end, a 
question remaining unsolved is what is the nature of the low-temperature peak in x'(T,u) 
for the D/J = 3.5 case, which will be cleared up as below. 

Figs. [3] and H] present the results of the temperature dependence of m z (T,uj), Xz(T,uj), 
and x'{T,u) from which one can see clearly evidences of MPT for the D/J = 3.5 case but 
not for the D/J = 10 case. m z (T, uj) is the averaged magnetization per spin projected along 
z direction, and Xz{T, uj) is the thermodynamic fluctuation of the magnetization. In Fig. [3l 
curves of m z (T, uj) are shown for these two cases with frequencies 1 x 10~ 4 < uj < 1 x 10~ 2 and 
H /J = 0.05. For the sake of reference, one curve of m z (T, uj) at uj = 1 x 10 -4 (i.e., the violet- 
colored solid line) is also plotted for D/J = 0, the case of non-anisotropic pure Heisenberg 
model possessing a well-known ferromagnetic phase transition [2lJ. For D/J = 3.5, the 
transition width of magnetization does not change until low frequencies uj < 5 x 10~ 4 with 
which the width gets narrower and narrower and m z (T,uj) curve approaches to the curve 
for D/J = 0. This change apparently corresponds to the change of the low-temperature 
peak with frequency in x'(T, ui) shown in Fig. [2j Strikingly, Xz{T,uj) in the inset of Fig. [3] 
exhibits a sharp peak similar to that of the D/J = CclSG, 1.6., db ferromagnetic-like MPT. 
This result supports the coexistence of MPT and GPT revealed for the case D / J = 4 of the 
RMA model in Eq. (CQ) [5j. In contrast, the transition in magnetization of the D/J = 10 
case is quite broad. This is indicated further in the inset of Fig. [3] by the noisy blurring peak 
of Xz{T, uj) whose height is orders of magnitude lower than those of D / J = 3.5 and D/J = 
cases. In addition, the magnetization at low temperatures for the D/J = 3.5 case is high in 
magnitude of 0.7, albeit smaller than 1.0 for the D/J = case, and frequency-independent 
against the small and chaotically frequency-dependent value of that for the D/J= 10 case. 
This feature is probably due to their different magnetic structures: the asperomagnet (known 
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in literature as a correlated spin-glass or a "ferromagnet" with wandering axis) in the former 
versus the speromagnet in the latter [HQ]- We believe that MPT for D/J = 10 likely does 
not exist or at least is smeared out by strong RMA and this is why the low-temperature 
peak is suppressed completely in x'{T, uj) curves. 

In general, a phase transition like those for the D/J = and D/J = 3.5 cases has been 
termed the dynamic transition, which is a true thermodynamic phase transition usually stud- 
ied together with the dynamic hysteresis in pure magnetic systems (the systems without any 



random defect or anisotropy to pin the magnetic domains) [23l . |24| . These phenomena occur 
due to a relaxational delay of the magnetization in response to the, say, oscillating field. 
When the oscillation period of the field is much less than the effective relaxation time of the 
magnetic system the hysteresis loop becomes asymmetric about the origin with a nonvan- 
ishing area and a " spontaneously broken symmetric phase" arises dynamically with a non- 
vanishing value of the dynamic order parameter Q, defined as Q(T, uj) = (u/2tt) § m(T, t)dt 
(Q(T,uj) is the period averaged magnetization and is equal to m z {T,uj) in our notation), 
where the instantaneous magnetization per site at time t and temperature T is calculated 
as m(T,t) = (l/L 3 )^2 i Si(T,t) ■ z. The system is in a dynamically-ordered phase when 
Q 7^ and the loop is asymmetric or in a dynamically-disordered phase when Q = and 
the loop is symmetric. A transition occurs at Td when one crosses the boundary separating 
the two phases. Notice that the boundary is dynamic in nature since Td depends on both 
H Q and uj, i.e., Td = Td(H ,uj). For any fixed frequency, the H Q — T plane is then divided 
by the dynamic phase boundary line Td(Ho,uj), which is in general convex towards the ori- 
gin. With large values of Hq, one gets a "forced oscillation" kind of scenario inducing the 
dynamically-disordered phase (Q = 0) at high T that gives rise to low values of T d (H ,uj). 
All of these features of the phase diagram have been obtained for pure magnetic systems 
of Ising models using mean-field and MC methods 24( . They are also observed in our MC 
simulation for the weak RMA Heisenberg model as shown in Fig. H]for the D/J = 3.5 case 
with Hq/ J = 0.05 ~ 2.5 and uj = 3 x 10~ 3 . Very interestingly, the dip in x'{T,u) as well as 
the only peak in x"{T, uj) occur somewhere around TAHq, uj) and quite similar to the fashion 



for two- and three-dimensional pure Ising models [24|. (note, however, that it is difficult 
to determine T d (H ,uj) in the m z (T,uj) (i.e., Q(T, uj)) curve because m z (T,uj) undergoes a 
gradually broad transition shown in Fig. HI Instead, we prefer to take the temperature T a 
at the peak in Xz{T,uj) or equivalently the temperature T p at the low-temperature peak in 
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x'(T, uj) to construct the diagram of the dynamic transition for the D/J = 3.5 case shown in 
the inset of Fig. [U) Eventually, in ac susceptibility measurements one may be indicated pre- 
cisely the same dynamic transition (where the peaks or dips are shown) as that the dynamic 
order parameter (if it could be directly measured in experiment) provides as long as the 
values of H are very small so that the dynamic transition is continuous because large values 
of Hq lead to a crossover of continuous/discontinuous transition at the tricritical point (not 
shown in our simulation) in the Hq — T diagram Therefore, the ac susceptibility for 

a-Ho28Fe72 [l^] is a particularly prominent example to study experimentally the dynamic 
transition in weak RMA systems using the ac susceptibility measurements. 

In summary, our MC simulation shows that the RMA model in Eq. ([T]) can be employed to 
understand the distinct behaviors in x'{T,cu) of a-Ho 2 8Fe 7 2 and Dy 40 Al24Co2oYiiZr 5 glasses 
Q, Q] where the nature of the low-temperature peak of x'{T, u) for the former is a dynamic 
transition. This result marks a striking similarity between weak RMA Heisenberg model 
and pure ferromagnetic spin models and sheds light on the nature of magnetic transition 
and magnetic ordering, i.e., QLRO, in magnets with weak random anisotropy. 

This work was financially supported by the National Science Council of Taiwan, R.O.C, 
under Grant No. NSC 97-2112-M-007-007-MY3. 
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1.5 




Temperature, 77 J 



FIG. 1: x'(T,oj) curves at u> = 3 x 10 3 and Hq = 0.05 for various D/J values. x'(T,u}) exhibits 
distinct characteristics of weak (D/J < 5) and strong (D/J > 5) degrees of RMA. 
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Temperature, T/J 



FIG. 2: Temperature dependence of x'(T,u) at various frequencies and Hq = 0.05 for D/J = 3.5 
(upper panel) and D/J = 10 (lower panel). Insets show laws of critical slowing-down dynamics for 
the two cases, respectively. 



10 




FIG. 3: Temperature dependence of m z (T,co) at various frequencies and Hq = 0.05 for D/J = 
(solid line), D/J = 3.5 (solid symbols), and D/J= 10 (open symbols). The inset shows Xz(T,o 
for these cases at w = 1 x 10 -4 . 
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0.8 




Temperature, T/J 



FIG. 4: Temperature dependence of m z (T,u), Xz(T,u), and x'(T,lu) for D/J = 3.5 at u = 3x 10 -3 
and different values of Ho (different colored curves). The one-peak and two-peak curves are Xz(T, uj) 
and x'(T,uj), respectively. The inset presents the dynamic phase diagram in Hq — T plane, where 
the T a line (or the T p line) shows the dynamic phase boundary between the dynamically-ordered 
phase (m z ^ 0) and dynamically-disordered phase (m z = 0). 



12 



